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Abstract We describe the complete mitochondrial genome of a small-scaled forest agamid (Pseudocalotes microlepis), 
which is 17 873 bp in size, containing 13 PCGs, 2 rRNAs, 22 tRNAs and non-coding regions. The mitogenome has a 
typical gene order among squamates. 13 PCGs include 2 start codons (ATG and ATA), 3 stop codons (TAG, TAA and 
AGG), and an incomplete stop codon (T-). Codon usage analyses showed that CUA-Leul and CGG-Arg are the most 
frequently and rarely used codon, respectively. All 22 tRNAs were predicted to form canonical cloverleaf secondary 
structures, except for two tRNAs (tRNA and tRNA" ^9?) lacking the dihydorouridine (DHU) arm. The large non- 
coding region (control region) is 2 687 bp long (28.3% T, 18.2% C, 42.3% A, and 11.2% G), with four different types 
of repeating sequences. The phylogenetic tree resulting from BEAST analyses based on concatenated 2 rRNAs and 13 
PCGs in sequence revealed that the newly sequenced P. microlepis, where the genus Acanthosaura were aggregated. 
Together with Calotes versicolor, they constitute the subfamily Draconinae. 
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1. Introduction 


The genus Pseudocalotes (Draconinae; Agamidae; 
Squamata) are Calotes-like lizards which inhabit 
mountain regions of Indo-China and the Sunda region, 
and are found mostly on trees or bushes in tropical 
mountain forests (Hallermann and Bóhme, 2000; Ziegler 
et al., 2006). To date, the complete mitochondrial 
genomes (mitogenomes) of 5.8% of Agamidae species 
(25/480) were available in GenBank, including 14 species 
in Agaminae, 4 species in Leiolepidinae, 3 species in 
Draconinae, 2 species in Amphibolurinae, 1 species in 
Hydrosaurinae and 1 species in Uromastycinae. However, 
the phylogenetic position and inter-relationships of the 
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subfamilies have yet to be determined. Some researchers 
proposed that the group (Agaminae + Draconinae) + 
Hydrosaurinae was the sister group of Amphibolurinae, 
and Uromastycinae was the outermost subfamily of 
agamids (Macey et al., 2000; Okajima and Kumazawa, 
2010). However, Pyron et al. (2013) proposed that the 
group of (Agaminae + Draconinae) was the sister group of 
(Amphibolurinae + Hydrosaurinae) with the study using 
5 nuclear loci (BDNF, c-mos, NT3, R35 and RAG-1) and 
5 mitochondrial loci (12S, 16S, Cytb, ND2 and ND4), 
and Uromastycinae was also the outermost subfamily 
of agamids, but the relationship Amphibolurinae + 
Hydrosaurinae was weakly supported. 

In this study, we sequenced the complete mitogenome 
of a small-scaled forest agamid (Pseudocalotes 
microlepis). This lizard occurs in the mountain forests of 
Hainan and Guizhou in China, Thailand, Laos, Myanmar 
and Vietnam (Ananjeva et al., 2011; Uetz, 2016; Zhao 
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and Adler, 1993). We analyzed the gene content, base 
composition, codon usage, tRNAs (transfer RNAs) 
structure and control region of this species. We then 
conducted a partitioned Bayesian phylogenetic analysis 
of related species based upon concatenated 2rRNAs 
(ribosomal RNAs) and 13 PCGs (protein-coding genes) 
in sequence (12S rRNA, 16S rRNA, ND1, ND2, COI, 
COL, ATP8, ATP6, CONI, ND3, ND4L, ND4, ND5, ND6 
and Cytb). We analysed the hitherto longest molecular 
data (14 024 bp) in Agamidae, and compared to the 25 
agamids with mitogenomes sequenced, in order to explore 
the phylogenetic relationships among the subfamilies. 


2. Materials and Methods 


2.1. Sample collection and DNA extraction The sample 
(voucher number XLHZ601) was collected from Hainan, 
China, and stored at —80°C in laboratory at Hangzhou 
Normal University. Total genomic DNA was extracted 
using the Genomic DNA kit (TransGen, China), according 
to the manufacturer-supplied protocols. 


2.2. Primer design, amplification and sequencing The 
species-specific primers were designed based on highly 
conserved sequences (Kumazawa and Endo, 2004), which 
were designed with software Primer Premier 5 and were 
identified using multiple alignments of the agamids (Table 
1). PCR was performed using a final reaction volume of 
25 uL, of 2.5 uL 10 x TransTaq HIFI Buffer, 2 uL dNTPs, 
1 uL forward and reverse primers for each, 0.5-2 uL 
DNA Template, 0.25-5 uL HIFI DNA Polymerase, and 
addition of double distilled water to a final volume of 25 
uL. The PCR procedure was conducted on a Mastercycler 
(Eppendorf, Germany) using the following program: pre- 
denaturation at 94°C for 5 min; 35 cycles of denaturation 
at 94°C for 30 s, annealing at 42°C—56°C for 30 s and 
extension at 72°C for 30 s; followed by a final elongation 
step of 72?C for 10 min. The PCR reaction products 
were electrophoresed in a 1.2% agarose gel and purified 
with PCR purification Kit (OMEGA, China). Then, each 
fragment of PCR was cloned into the pEASY-T5 Zero 
Cloning vector (TransGen, China) and sequenced with 
M13 primers from both directions by a primer-walking 
strategy. Each sequence overlapped the next contig 
between 150—300 bp. 


2.3. Sequence analysis We analyzed DNA sequences 
and performed contig assembly with the software Seqman 
(DNASTAR). We identified 13 PCGs using ORF Finder 
implemented at the NCBI website (https://www.ncbi.nlm. 
nih.gov/orffinder/). We used tRNAscan-SE search server 
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(Lowe and Eddy, 1997) and MITOS web servers (Bernt et 
al., 2013) to predict the secondary structure and anticodon 
sequences of tRNAs. We determined the boundaries 
of 2 rRNAs based on alignments of rRNA sequences 
of other species, as with Calotes versicolor (Amer and 
Kumazawa, 2007). We then compared our sequences 
with the available agamid mitogenomes using Clustal 
X 2.0 (Larkin et al., 2007; Thompson et al., 1997), and 
searched the tandem repeat sequences of control region 
using Tandem Repeats Finder 4.0 (Benson, 1999). Then 
we calculated the GC and AT skew respectively following 
Perna and Kocher (1995) formula: AT skew = (A-T) / 
(A+T) and GC skew = (G-C) / (G+C) (Perna and Kocher, 
1995). We analyzed the nucleotide composition and 
relative synonymous codon usage (RSCU) values with 
Mega 6.0 (Tamura et al., 2013). The complete mtDNA 
sequence of P. microlepis was deposited to GenBank 
under accession number (KX898 132). 


2.4. Phylogenetic analyses To construct a phylogenetic 
tree, we used all 25 available mitogenomes in Agamidae, 
and used Chamaeleo calyptratus (Chamaeleonidae) as 
the outgroup. Accession number and the whole size for 
all mitogenomes were presented in Appendix 1. We used 
the nucleotide sequences consisted of 15 mitochondrial 
genes, and deleted the start and end codons of PCGs 
from analyses. Phylogenetic analyses were conducted 
using Bayesian uncorrelated lognormal approach in 
BEAST 2.0 (Bouckaert et al., 2014). Prior to estimating 
BEAST, we carried out the best partition schemes and 
corresponding nucleotide substitution models for each 
partition, using PartitionFinder 1.1.1 (Lanfear et al., 
2012). The best-fitting model was determined using the 
Bayesian Information Criterion (BIC). The partitions 
and models were listed in Table 2. We used one fossil 
calibration point, the estimated age of the split between 
oviparous and viviparous species of Phrynocephalus: 
9.73 (95% interval: 7.21—13.04) Ma (million years ago) 
(Jin and Brown, 2013). We used the relaxed lognormal 
clock model, and specified the standard Yule speciation 
process for the tree. Two independent runs of four heated 
MCMC chains (three hot chains and one cold chain) 
were simultaneously run for 200 million generations, 
with sampling conducted every 10,000 generations 
(Fitze et al., 2011; Jin et al., 2015; Kyriazi et al., 2008). 
We compared the results of two independent runs with 
Tracer 2.2.1 (Rambaut and Drummond, 2007). Then we 
discarded the first one million trees as “burn-in”. All four 
chains achieved the recommended adequate effective 
sample size of 200 for likelihood (Drummond et al., 
2006; Lin and Wiens, 2017). 
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Table 1 Primers designed to amplify and sequence in this study. L and H refer to forward primer and reverse primer, respectively; Tm refers 
to annealing temperature. 


Primer name Primer sequence (5'-3') Tm (*C) Reference 
112S-1L AGGATTAGATACCCTACTA 48 Kumazawa and Endo, 2004 
rl6S-3H CAKKTTTCCCTTGCGGTACT 

ul2S-1L GCGYACAYAY CGCCCGTC 48 Kumazawa and Endo, 2004 
16sbr-H CCGGTCTGAACTCAGATCACGT 

16SL CCAACAAGCCAACCGACA 52 This study 

16SH CGCCCATGAGTGCGTATT 

16S-NDIL TACGACCTCGATGTTGGATCAGG 46 This study 
16S-NDIH ATTTTTCGTAGTTGGGTTTGRTT 

NDIL CTAACAATAAAAGCGGGCAC 52 This study 

NDIH GGCTAGATGTAGGGCAAAGAT 

NDI-COIL GCCTGATACACTGGGACCT 54 This study 
ND1-COIH TGTTCATCCTGTGCCGAC 

rCOI-IL ATCGGCGGRTTYGGAAACTG 50 Kumazawa and Endo, 2004 
rCO2-1H TGGAAGTGWARTAGYTCTTCTAT 

COIIL CCMCGACTACTAGARGTAGA 53 This study 

COIIH GYGAATACGTAGGCTTGRAT 

uLys-1L AGCACTAGCCTTTTAAGC 46 This study 
uCO3-1H AAYGTCTCGTCATCATTG 

CO3-ND3L ATAGTWGACCCMAGCCCATGACC 56 Kumazawa and Endo, 2004 
CO3-ND3H GGGTCRAAKCCRCATTCRTA 

rND4-3L CCAAAAGCCCAYGTAGARGC 54 Kumazawa and Endo, 2004 
rCUN-IH CTTTTACTTGGADTTGCACC 

ND3-NDSL CGCACTCCAGGCAATAGAA 53 This study 
ND3-ND5H ACGGCAAGTTTGGGGTTA 

NDSL CCTGACACCCGAAAACGA 56 This study 

ND5H TGTGATGGAGGCTAAATGAAGT 

2NDSL TGGGGTCACTTCCACTCTTT 56 This study 
2ND5H GCGTTTGATGTGAGGGGTAT 

rHis-2L AACAAAAACAYTAGRCTGTG 42 Kumazawa and Endo, 2004 
rND5-2H ATWGYGTCTTTTGAGTARAAKCC 

rND6-3L GCAACWGAATAHGCAAATAC 44 Kumazawa and Endo, 2004 
rcytb-1H GCGTAGGCRAATAGGAAGTATCA 

ND6-CBL GARCARGACATYCGAAAAATRGG 50 This study 
ND6-CBH CAGCCGTGTTGTACATCTCG 

CBL TTGCCCTTTCAACCTCT 50 This study 

CBH GCGGCTCTTATTTGTATTC 

CB2L AGCCAACTGGGMCTMATAAT 42 This study 

CB2H GGTTTACAAGACCAGTTGCTTT 

CR2L CAAACCAGAGTGATACTTTC 50 This study 

CR2H CCCTCTAGGTTGATGTAGA 

CR3L CTTACATGAATCGGAAGC 48 This study 

CR3H GGCATATAGCCTCTTGAT 


Table 2 The best partition schemes and nucleotide substitution models for mitochondrial data based on BEAST carried out using 
partitionFinder. 


Subset Best Model Subset Partitions 

Pl GTR+I+G 12S, 16S, ATP6_pos1 

P2 GTR+HG COIL posl, NDI posl, ND3 posl, Cytb posl 

P3 GTR+HG COIL pos2, COII pos2, COI pos2, ND1_pos2 

P4 GTR+G ATP6 pos3, ATP8 pos3, ND1_pos3, ND2 pos3, ND3 pos3, ND4L_pos3, ND4 pos3, ND5 pos3, Cytb_pos3 
P5 GTR+HG ATPS8 posl, ND2 posl, NDAL posl, ND4 posl, ND5 posl 

P6 GTR+G ATP6 pos2, ATP8 pos2, ND2_pos2, ND3 pos2, NDAL pos2, ND4 pos2, ND5 pos2, Cytb pos2 

P7 SYM+I+G COII posl, COL posl 

P8 HKY+G COII pos3, COI pos3, COL pos3 

P9 TrN-I-G ND6 posl, ND6 pos2 


P10 HKY+G ND6_pos3 
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In order to estimate the substitution rate of each gene, 
along with their confidence intervals, we performed an 
additional BEAST analysis with clock models linked 
by 'gene' and nucleotide substitution models unlinked. 
Models of nucleotide evolution of each gene partition 
were calculated in jModelTest 2.1.7 (Darriba et al., 
2012), under the Akaike Information Criterion (AICc). 
Phylograms were drawn with Figtree 1.4. 


3. Results and Discussion 


3.1. Genome organization and structure The 
mitogenome of P. microlepis 1s a typical circular DNA 
molecule of 17 873 bp in length, similar in size to the 
other available mitogenomes of species in the Agamidae. 
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In comparison with the other agamids, the mitogenome 
of P. microlepis is longer than all species (Appendix 1) 
except Phrynocealus axillaris (17 937 bp). The difference 
in size is mainly due to the variable number of tandem 
repeats (VNTRs) in the control region. It contains a 
typical set of gene content: 13 PCGs, 2 rRNAs, 22 tRNAs 
and non-coding regions. Among these, 29 genes (12 
PCGs, 15 tRNAs and two rRNAs) are located on heavy 
(H) strand, and other genes (ND6 and seven tRNAs) are 
located on light (L) strand (Table 3). Gene overlaps of 
42 bp have been found at 10 gene junctions, the longest 
overlap (10 bp) exists between ATP8 and ATP6. The gene 
order of the P. microlepis mitogenome is identical to that 
of most squamates (Macey et al., 2006; Ujvari et al., 
2007). 


Table 3 Mitochondrial genome organization of Pseudocalotes microlepis. L and H refer to forward primer and reverse primer, respectively. 


Gene Strand Position Size(bp) Start codon Stop codon  Intergenic spacer 
tRNA™ H 1-67 67 0 
12S rRNA H 68-913 846 0 
tRNA” H 914-979 66 -1 
16S rRNA H 979-2 497 1519 1 
tRNA "Y H 2 499-2 570 72 2 
NDI H 2 573-3 541 969 ATG TAG E 
tRNA?" L 3 538-3 609 72 5 
tRNA” H 3 615-3 681 67 -l 
tRNA™* H 3 681-3 747 67 0 
ND2 H 3 748-4 776 1029 ATG TAA 5 
tRNA’? H 4 782-4 854 73 2 
tRNA“" L 4 857-4 924 68 4 
tRNA^" L 4 929—5 000 72 23 
tRNA“ E 5 024-5 076 53 -1 
tRNA™ L 5 076-5 139 64 0 
COI H 5 140—6 687 1548 ATG TAG 29 
tRNA?" UC» È 6 717-6 786 70 4 
tRNA“? H 6 791-6 858 68 3 
COII H 6 862-7 548 687 ATG AGG -7 
tRNA’ H 7 542-7 609 68 1 
ATP8 H 7 611-7772 162 ATG TAA -10 
ATP6 H 7 763-8 444 682 ATG T 0 
COIII H 8 445-9 228 784 ATG T 0 
tRNA?* H 9 229-9 294 66 2 
ND3 H 9 297-9 641 345 ATG TAA -2 
tRNA“® H 9 640-9 706 67 1 
ND4L H 9 708-9 998 291 ATG TAA -7 
ND4 H 9 992-11 359 1368 ATG AGA 7 
tRNA" H 11 367-11 431 65 

tRNA9" ^69 H 11 432-11 489 58 8 
tRNA H 11 498-11 569 72 

ND5 H 11 570-13 348 1779 ATA TAG -8 
ND6 L 13 341-13 853 513 ATG AGG 0 
tRNA“ L 13 854-13 920 67 2 
Cytb H 13 923-15 048 1126 ATG T 9 
tRNA™ H 15 058-15 122 65 “i 
tRNA’? H 15 122-15 186 65 0 
control region H 15 187-17 873 2687 
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3.2. Nucleotide composition Similar to most other 
mitogenomes in Agamidae, the nucleotide composition 
of P. microlepis mtDNA is biased toward A and T. The 
overall A + T content of mitogenome is 59.2% (35.3% 
A, 23.9% T, 27.6% C and 13.2% G). The AT skew and 
GC skew is 0.1943 and —0.3541, respectively. The P. 
microlepis mitogenome has a distinct bias against G at 
first codon position (A: 37.0%, T: 22.5%, C: 27.2% and 
G: 13.3%). The percentage of purines (48.8%) is slightly 
lower than pyrimidines (51.2%) at the second position 
and the third position. 


3.3. Protein-coding genes and relative synonymous 
codon usage The total length of the 13 PCGs in 
P. microlepis mitogenome is 11 283 bp, accounting for 
63.13% of the entire mitogenome sequence. All the PCGs 
initiated with a typical start codon (ATG), except ND5 
which starts with ATA. Among stop codons, TAA is the 
most common. ND2, ATP8, ND3 and ND4L end with 
TAA; NDI, COI and NDS end with TAG; COII and ND6 
end with AGG; ND4 ends with AGA; ATP6, COIII and 
Cytb end with an incomplete end codon (T-). The post- 
transcriptional polyadenylation can produce a standard 
TAA stop codon (Han and Zhou, 2005). 

The relative synonymous codon usage (RSCU) for 
the mitochondrial PCGs in P. microlepis exhibited 62 
amino-acid encoding codons as well as an over-usage of 
A and T at the third codon positions (Table 4). Among 
them, CUA-Leul (7.12), ACA-Thr (6.02) and AUA-Met 
(5.19) are the most frequently used codons. The least- 
frequent codons are CGG-Arg (0.03), CGU-Arg (0.05) 
and UCG-Ser2 (0.13). These codons are composed of A 
and U nucleotides, indicating a high usage of A and T in 
P. microlepis PCGs. 


3.4. Ribosomal and transfer RNA genes The 2 rRNAs 
(12S and 16S) of P. microlepis are located between 
tRNA™ and tRNA‘ 9 and separated by tRNA"' 
(Table 3). The lengths of 12S rRNA and 16S rRNA are 
determined to be 846 bp and 1 519 bp, and it varies from 
830 bp in A. lepidogaster to 930 bp in A. armata. 16S 
rRNA varies from 1 479 bp in genus Phrynocealus to 
1 567 bp in Leiolepis boehmei, respectively. The size 
is similar to that of other metazoan mtDNA (Zhang et 
al., 2009). The typical set of 22 tRNA of P. microlepis 
is ranging in size from 53 bp for tRNA“ to 73 bp for 
tRNA", as similar to other metazoan mitogenomes 
(Yoon et al., 2015). The 22 tRNAs possess a canonical 
cloverleaf secondary structure composed of four arms 
(dihydorouridine arm, anticodon arm, TC arm and 
aminoacyl acceptor arm) with conserved size (Figure 1). 


Table 4 Codon usage in Pseudocalotes microlepis mitochondrial 
protein-coding genes. A total of 3 737 codons for analyzed, 
excluding the start and stop codons. AA, amino acid; RSCU, 
relative synonymous codon usage; n = frequency of each codon; 
% = n/3737. 


AA Codon n % RSCU 
Phe (F) UUU 78 2.09 0.83 
UUC 109 2.92 LiT 
Leu2 (L2) UUA 87 2.33 0.88 
UUG 24 0.64 0.24 
Leul (L1) CUU 97 2.6 0.98 
CUC 80 2.14 0.81 
CUA 266 7.12 2.7 
CUG 37 0.99 0.38 
Ile (T) AUU 125 3.34 0.9 
AUC 154 4.12 L1 
Met (M) AUA 194 5.19 1.62 
AUG 46 1.23 0.38 
Val (V) GUU 37 0.99 0.93 
GUC 24 0.64 0.6 
GUA 73 .95 1.84 
GUG 25 0.67 0.63 
Ser2 (S2) UCU 43 .15 0.89 
UCC 73 95 1.51 
UCA 120 321 2.47 
UCG 5 0.13 0.1 
Pro (P) CGU 15 0.4 0.3 
CCC 38 02 0.76 
CCA 131 3.51 2.61 
CCG 17 0.45 0.34 
Thr (T) ACU 40 .07 0.36 
ACC 161 4.31 1.46 
ACA 225 6.02 2.05 
ACG 14 0.37 0.13 
Ala (A) GCU 36 0.96 0.55 
GCC 129 3.45 1.95 
GCA 86 2.3 1.3 
GCG 13 0.35 0.2 
Tyr (Y) UAU 26 0.7 0.5 
UAC 79 2.11 1.5 
Stop UAA 0 0 0 
UAG 0 0 0 
His (H) CAU 14 0.37 0.3 
CAC 80 2.14 1.7 
Gln (Q) CAA 84 2.25 1.7 
CAG 15 0.4 0.3 
Asn (N) AAU 30 0.8 0.39 
AAC 122 3.26 1.61 
Lys (K) AAA 111 2.97 1.87 
AAG 8 0.21 0.13 
Asp (D) GAU 15 0.4 0.55 
GAC 40 1.07 1.45 
Glu (E) GAA 58 1.55 1.47 
GAG 21 0.56 0.53 
Cys (C) UGU 9 0.24 0.72 
UGC 16 0.43 1.28 
Trp (W) UGA 83 221 1.75 
UGG 12 0.32 0.25 
Arg (R) CGU 2 0.05 0.13 
CGC 9 0.24 0.59 
CGA 49 1.31 321 
CGG 1 0.03 0.07 
Ser (S1) AGU 16 0.43 0.33 
AGC 34 0.91 0.7 
Stop AGA 0 0 0 
AGG 0 0 0 
Gly (G) GGU 20 0.54 0.4 
GGC 54 1.45 1.07 
GGA 82 2.19 1.63 
GGG 45 1.2 0.9 
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c AA A C  AAAC A e AAGC c , AAAC A UUGC U 
UA A A A A Cy A A A c G GA G G 
U-AA U-AA U-AA C-GC G AG 
A-U U-A U-A U-A U-A 
G-c A-U A-U eeg U-A 
c-G G-C G-C Geg G-C 
A-U G-C A A e£ G. U 
U A c c c A U U U G 
U A c A U A c A U A 
U G 
GA Gys Ug U G UGA 
Leucine y : Valine EN 
(Leu, L2) C-G (Val, V) A-u 
(CUN) U-À A-U 3 
U-A A-v———Aminoacyl acceptor arm 
V n 
= A-U 
AP "cance dus FER 
CA 
A ga M, e i 
: ERI i Uc GUUAA , 
M T TYC arm 


pd — —— Variable loop 
G-C 
CHG Anticodon arm 
c A 
U G 
U A G 


Dihydorouridine arm 


Figure 1 Putative 22 tRNAs secondary structures of Pseudocalotes microlepis. The minus (-) indicates Watson-Crick base pairing, and dots 
indicate G-U base pairing. It is composed of Aminoacyl acceptor (AA) arm, Dihydorouridine (DHU) arm, Anticodon (AC) arm, TYC (T) 
arm and Variable loop. 
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Whereas two tRNAs (tRNA“ and tRNA" ^9?) appear 
to lack the dihydorouridine (DHU) arm. The loss of the 
DHU arm in tRNA?" ^9? has been considered a common 
condition of metazoan mitogenomes (Wolstenholme, 
1992). However, the loss of the DHU arm in tRNA“ is 
an unusual phenomenon, which has also been observed 
in Gekko gecko (Han and Zhou, 2005). Further research 
is needed to determine the molecular mechanisms 
responsible for keeping such defective tRNAs functional. 


3.5. Non-coding regions The small non-coding region 
includes several intergenic spacers, ranging from 1 
to 29 bp (Table 3), most of which are shorter than ten 
nucleotides. The longest intergenic spacer sequence we 
found is located between COI and tRNA?*" '*9, The 
large non-coding region (control region) is 2 687 bp in 
size and located between tRNA” and tRNA”. The size 
is remarkably longer than other species in Draconinae 
because of the VNTRs. The nucleotide composition 
is 42.3% A, 28.3% T, 18.2% C and 11.2% G, with a 
strong bias use of G. The structure is typical including 
Termination-Associated Sequence (TAS) and Conserved 
Sequence Blocks (CSB) (Jin et al., 2015; Shi et al., 2013; 
Xiong et al., 2010). VNTRs contain four distinct tandem 
repeat units (15 950-16 014, 16 018—16 086, 16 138-16 
671 and 16 707—17 857). They are 65 bp, 43 bp, 534 bp, 
and 1151 bp in length, respectively. The small tandem 


70 60 50 40 30 20 
Million years ago 


repeats units as 5'-AACA-3' and 5'-A/G (G) CAA-3' 
have 16.3 copies and 10.8 copies, respectively. One large 
tandem repeats (74 bp) have 7.2 copies, another one (75 
bp) have 15.4 copies. VNTRs have also been regarded as 
a common feature for the mitogenomes of reptiles (Xu 
and Fang, 2006), and could provide reliable phylogenetic 
and genetic information for closely related species 
(Zardoya and Meyer, 1998). 


3.6. Phylogenetic analysis 
relationships were constructed using BEAST based on 


The phylogenetic 


15 mitochondrial genes. The final alignment resulted in 
14 024 nucleotide sites for 26 ingroup and one outgroup 
taxa. The number of sequences and substitution rates, 
multiple sequence alignments length, and models of genes 
were reported for each gene in Table 5. The topology of 
phylogenetic tree was shown in Figure 2. Most nodes 
were well supported by high posterior probabilities. The 
divergence between Chamaeleonidae and Agamidae 
was estimated at 64.87 Ma; within Agamidae, the 
basal branching split was estimated at 60.02 Ma. The 
divergence between oviparous and viviparous species of 
Phrynocephalus was 9.20 Ma. 

Our results revealed that the newly sequenced 
P. microlepis and the genus Acanthosaura were 
aggregated, and together with C. versicolor they 
constitute the subfamily Draconinae. However, the 


Chamaeleo calyptratus 
Pseudotrapelus sinaitus 
Xenagama taylori 
Ir— Phrynocephalus albolineatus 
1 Phrynocephalus grumgrzimailoi 
Phrynocephalus przewalskii 
1 Phrynocephalus versicolor 
Phrynocephalus helioscopus 
Phrynocephalus axillaris 
1lo.56 Phrynocephalus mystaceus 
0.96 Phrynocephalus erythrurus 
l Phrynocephalus theobaldi 
1 lf Phrynocephalus guinanensis 
Phrynocephalus putjatai 


Phrynocephalus forsythii 
1 Acánthosaura armata 


Acanthosaura lepidogaster 
Pseudocalotes microlepis 
Calotes versicolor 
Hydrosaurus amboinensis | Hydrosaurinae 
Chlamydosaurus kingii ; ; 
Pogona vitticeps id | Amphibolurinae 
l Leiolepis belliana 
Leiolepis reevesii 
Ir Leiolepis boehmei 
Leiolepis guttata 
Uromastyx benti | Uromastycinae 


Agaminae 


Draconinae 


Leiolepidinae 


Figure 2 Time calibrated Bayesian Phylogenetic tree of amagids and one outgroup estimated using BEAST based on mitochondrial genes 
(concatenated 2 rRNAs and 13 PCGs) for Markov chains. Numbers on nodes are posterior probability values. Hatched rectangles indicate 


95% credibility range for divergence times. 
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Table 5 Number of sequences and substitution rates, length of the gene fragments, models of genes for each gene as selected by jModelTest 


according to the AICc. 


Substitution rates 


Gene Length (bp) Best model -InL (% of substitutions per site per million years per lineage) 

Mean Stderr of mean 9595 HPD Interval 
128 1048 GTR+G 11069.069 0.7536 4.50E-4 (0.6791, 0.8315) 
16S 1612 GTR+HG 15784.0224 0.5624 3.15E-4 (0.5132, 0.6085) 
NDI 948 GTR+I+G 11729.8923 1.0424 6.77E-4 (0.9283, 1.1605) 
ND2 1032 GTR+I+G 14358.5244 1.1804 6.60E-4 (1.0703, 1,2942) 
COI 1539 GTR+I+G 16784.6191 0.9977 6.67E-4 (0.8984, 1.0992) 
COII 672 GTR+I+G 7982.3681 0.9795 8.78E-4 (0.8484, 1.1165) 
ATP8 156 HKY+G 2460.6986 1.0932 1.25E-3 (0.8958, 1.3291) 
ATP6 687 GTR+I+G 9884.5279 1.2363 8.39E-4 (1.0947, 1.3743) 
COIII 780 GTR+I+G 8866.9065 1.0308 8.27E-4 (0.8964, 1.1679) 
ND3 339 GTR+I+G 4535.8185 1.0638 1.27E-3 (0.8916, 1.2598) 
ND4L 291 HKY+I+G 4155.6149 1.006 8.95E-4 (0.8603, 1.1685) 
ND4 1365 GTRH+G 18320.3413 1.0761 5.67E-4 (0.9871, 1.1636) 
ND5 1845 GTRH+G 26176.2835 1.134 4.61E-4 (1.0596, 1.2062) 
ND6 567 GTRH+G 7683.8142 1.1316 7.22E-4 (1.0112, 1.2584) 
Cytb 1143 GTR+I+G 15413.9802 1.0915 6.30E-4 (0.9867, 1.1908) 
usage of mitogeneome did not allow us to resolve with References 


support the position of Hydrosaurinae which was instable 
across previous studies. For example, some previous 
studies (Blankers et al., 2013; Townsend et al., 2011; 
Wiens et al., 2012) placed Hydrosaurinae as sister to 
Amphibolurinae + (Agaminae + Draconinae), other 
studies (Macey et al., 2000; Okajima and Kumazawa, 
2010) placed Hydrosaurinae as sister to (Agaminae + 
Draconinae), whereas we and Pyron et al. (2013) placed 
Hydrosaurinae as the sister-group to Amphibolurinae with 
weak support. Further studies were required to resolve the 
position of Hydrosaurinae. 


4. Conclusions 


In this study, we sequenced and annotated the complete 
mitogenome of P. microlepis. Our results present the gene 
content, base composition, codon usage, tRNAs structure, 
VNTRs in the control region and phylogenetic analysis of 
related species. This is the first complete mitogonome of 
the genus Pseudocalotes. The research is intended to be 
helpful for the exploration on the phylogenetic position 
and interrelationships of the subfamilies in Agamidae. 
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Appendix 1 Species used in phylogenetic ananlysis. 


(Sub)family 
Agamidae 
Agaminae 


Amphibolurinae 


Draconinae 


Hydrosaurinae 
Leiolepidinae 


Uromastycinae 
Chamaeleonidae 


Species 


Phrynocephalus albolineatus 
Phrynocephalus axillaris 
Phrynocephalus erythrurus 
Phrynocephalus forsythii 
Phrynocephalus grumgrzimailoi 
Phrynocephalus guinanensis 
Phrynocephalus helioscopus 
Phrynocephalus mystaceus 
Phrynocephalus przewalskii 
Phrynocephalus putjatai 
Phrynocephalus theobaldi 
Phrynocephalus versicolor 
Pseudotrapelus sinaitus 
Xenagama taylori 
Chlamydosaurus kingii 
Pogona vitticeps 
Acanthosaura armata 
Acanthosaura lepidogaster 
Calotes versicolor 
Pseudocalotes microlepis 
Hydrosaurus amboinensis 
Leiolepis belliana 
Leiolepis boehmei 
Leiolepis guttata 
Leiolepis reevesii 
Uromastyx benti 
Chamaeleo calyptratus 


Size (bp) 


16 808 
17 937 
16 379 
16 143 
16 301 
16 279 
16 249 
16 660 
16 892 
16 283 
16 608 
16 429 
16 560 
16 220 
16 761 
16 751 
16 544 
16 899 
16 670 
17 873 
16 129 
16 109 
16 106 
16 552 
16 908 
16 380 
17 433 


Accession No. 


KP279760 
KC119493 
KJ630904 
KP232959 

KM093859 
KJ885621 

KM093858 
KC578685 
KF572032 
KJ830752 
KJ551842 
KJ749841 
AB262447 
DQ008215 
EF090421 
AB166795 
AB266452 
KR092427 
AB183287 
KX898132 
AB475096 
AB537554 
AB537555 
AB476400 
KJ530718 
AB114447 
EF222192 


Source 


Shao et al., 2015 
Li et al., 2013 
Zhu et al., 2016 
Shao et al., 2016 
Shuang et al., 2016 
Fu et al., 2016 
Li et al., 2016 
Chen et al., 2014 
Li et al., 2015 
Tong and Jin, 2016 
Liao and Jin, 2016 
Song et al., 2016 
Okajima and Kumazawa, 2010 
Macey et al., 2006 
Ujvari et al., 2007 
Amer and Kumazawa, 2005 
Okajima and Kumazawa, 2010 
Yu et al., 2015 
Amer and Kumazawa, 2007 
This study 
Okajima and Kumazawa, 2010 
Srikulnath et al., 2009 
Srikulnath et al., 2009 
Okajima and Kumazawa, 2010 
Tong et al., 2016 
Okajima and Kumazawa, 2010 
Macey et al., 2008 


